Brownian ratchets driven by asymmetric nucleation of hydrolysis waves 
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We propose a stochastic process wherein molecular transport is mediated by asymmetric nucle- 
ation of domains on a one-dimensional substrate. Track-driven mechanisms of molecular transport 
arise in biophysical applications such as HoUiday junction positioning and coUagenase processivity. 
In contrast to molecular motors that hydrolyze nucleotide triphosphates and undergo a local molec- 
ular conformational change, we show that asymmetric nucleation of hydrolysis waves on a track can 
also result in directed motion of an attached particle. Asymmetrically cooperative kinetics between 
"hydrolyzed" and "unhydrolyzed" states on each lattice site generate moving domain walls that 
push a particle sitting on the track. We use a novel fluctuating- frame, finite-segment mean field 
theory to accurately compute steady-state velocities of the driven particle and to discover parameter 
regimes which yield maximal domain wall flux, leading to optimal particle drift. 

PACS numbers: 82.39.-k,87.16.Ac,05.40.-a 
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Molecular motors such as kincsins, myosins, helicascs, 
and polymerases typically convert part of the free energy 
of ATP or GTP hydrolysis to a conformational change 
This molecular deformation leads to motion of the 
motor against a load on a track. Although the literature 
on such molecular motors is vast, much less attention 
has been paid to the theory of molecular motions that 
exploit the dynamics of the track on which translation 
occurs. Such loads are propelled by the track, which itself 
is undergoing catalyzed state changes by, e.g., hydrolysis. 

Two biological strategies involving track-propelled par- 
ticles are collagenase catalysis and Holliday junction 
transport. Collagenase MMP-1, an enzyme that asso- 
ciates with and cleaves collagen, is propelled by proteol- 
ysis of the collagen track 0, 0] . The cleaving of bonds 
prevents the collagenase from diffusing back across the 
broken bond, resulting in biased transport of the collage- 
nase. Thus the statistical dynamics of the track propels 
the enzyme. This dynamic has been modeled by a burnt 
bridge model [1,0,11,11. 

Another system where substrate modification possibly 
leads to biased motion is the translocation of Holliday 
junctions 0. The Holliday junction at which two double- 
stranded DNA molecules exchange one of their strands 
may be moved by the dynamics of hydrolysis states of 
the DNA binding protein RecA. RecA polymerizes on at 
least one of the strand-exchanging dsDNA molecules, as- 
sembling into a long nucleoprotein filament. The RecA 
monomers appear to hydrolyze ATP and can exist in dif- 
ferent states, much like the intermediate hydrolysis states 
of myosin motors. The dynamics of the interconversion 
among these hydrolysis states may provide the force nec- 
essary to rotate DNA strands about each other during 
Holliday junction translocation. An especially promising 
model of this process exploits asymmetric cooperativity 
in the hydrolysis of the nucleotide triphosphate cofac- 
tors associated with each RecA monomer. This intrinsic 
asymmetry of the filament gives rise to "waves" of hy- 



drolyzed monomers with a preferred direction, thereby 
moving the junction by virtue of its preferential attach- 
ment to the hydrolyzed segment of the RecA filament 
01 . These examples constitute only two of many mecha- 
nisms through which chemical energy may be harnessed 
to perform mechanical work by the substrate rather than 
a motor protein. In this Letter, we develop a gen- 
eral stochastic theory of track-driven, hydrolysis wave- 
mediated transport. In addition to analyzing our model 
using Monte-Carlo (MC) simulations, we also formulate 
a moving-frame mean field theory (MFT) that accurately 
predicts novel features of the transport. 
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FIG. 1: Schematic of the asymmetric nucleation process. An 
intrinsic asymmetry in the lattice sites gives rise to asymmet- 
ric cooperativity and nucleation. The transported particle is 
represented by a triangle. 

As in models of the ATP cycle of myosin motors, in 
RecA hydrolysis wave-mediated transport the RecA sub- 
units can exist in a number of substates corresponding 
to sites that have bound ATP, ADP-|-Pi, ADP, or are 
empty. To simplify our model, we will assume that each 
site of the substrate lattice exists in only one of two pos- 
sible states, "hydrolyzed" (cr = 1) and "unhydrolyzed" 
(a = 0). Any lattice site i can transition from state 
CTi = 1 to state (Ti = with rate kg. The reverse process, 
physically corresponding to "hydrolysis" or "nucleation," 
fills an empty site. In our model, an asymmetry arises 
in the nucleation transitions ai = crt ~ 1. If site 
I — 1 is also in the state Ci-i = 0, then the transition 
(T; = ^ (Ti = 1 occurs with rate fc_. However, if 
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ai-i = 1, then the transition cr^ = — > di = 1 occurs 
with rate . If ^ k- , the process is asymmetric and 
can lead to a steady-state current of domain walls. If a 
particle is associated with the lattice, it will be pushed 
each time a domain wall passes it. Thus, a net flux of do- 
main walls will lead to directed particle transport. The 
kinetics of the lattice is outlined in Fig. [1] 

The corresponding Master equation is similar to that 
which describes Glauber dynamics of a one-dimensional 
Ising model, except that the asymmetry in the transi- 
tion rates prevents this system from supporting an equi- 
librium state. Since no exact solutions are known, wc 
employ a hybrid finite-segment mean field theory (MFT) 
and MC simulations to obtain quantitative results and 
physical understanding. 

First consider a translationally invariant (infinite or 
periodic) lattice in the absence of an associated load 
particle. In the fixed laboratory frame, the moments 
{(TiCj . . .) of the hydrolysis states can be derived from 
the Master equation. These hierarchical moment equa- 
tions are not closed. For example, the equation for the 
first moment 



dt 



k- - (fc_ + /co)((7,) + A(<T,_i) - A(cricr,_i), 



(1) 

where A = k^ — fc_ is the hydrolysis asymmetry, depends 
on correlations {aiUi-i). The simplest mean field ap- 
proximation assumes {uiGi-i) = {ai-i){ai), which, when 
combined with the steady-state limit {A{(Ti)/dt = 0) of 
Eq. [1] gives a steady-state mean hydrolysis level 



2fc_ - fco A/(fc7^^oF+40?o 



(2) 
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As shown by the dotted curves in Fig. [2^, this result is 
only in qualitative agreement with the mean hydrolysis 
level obtained from MC simulations on a lattice with N = 
1000 sites(open circles). 

The locality of the asymmetric interactions suggests 
that correlations are short-ranged as in the totally asym- 
metric exclusion process Q. Thus, more accurate 
approximations can be systematically implemented by 
considering small clusters in which all possible config- 
urations are identified, and enumerating all the transi- 
tions among them. The densities at both ends of this 
cluster are then self-consistently matched to bulk values 
inferred by the statistics within the cluster. This finite- 
segment mean field approach has been used to study the 
nonequilibrium steady-states of related models such as 
the asymmetric exclusion process [13, U | . For example, 
consider all possible configurations in a segment of m = 2 
lattice sites. If we enumerate the states corresponding 
to the binary representation of the state number, {i.e 



Po = 00, Pi = 01, P2 = 10, P3 = 11) 
tion matrix defined by P = MP is 



the 2^ X 2^ transi- 
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(3) 

where P = (Pq, Pi, P2, Ps)"^ is the probability vector and 
s is the mean occupancy of the site immediately to the 
left of the explicitly enumerated pair. Since s repre- 
sents the mean occupation of the rightmost site of the 
preceding segment, we impose self-consistency by setting 
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Pi = s and solving for s numerically. The simple 



mean field approximation (Eq. ^ corresponds to m = 1. 
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FIG. 2: (a) Mean densities computed from finite-segment 
MFT and from MC simulation. The three groups of den- 
sities correspond to k- — 0.1, 1, 10. Within each group, 
finite-segment mean field results with sizes m = 1, 2, 3, 4 (de- 
noted MFTm) are compared with results from MC simula- 
tions (open circles), (b) Density profiles (for k+ = 1.8,3) in 
the lab frame (thin light curves) and in a frame moving with 
the load particle (thick dark curves) that follows the rules 
indicated in Fig. [S] 

Fig. shows the increasing accuracy in determin- 
ing (a) upon using larger m (shown are m = 1,2,3,4) 
in the finite-segment mean field approach. Although 
m = 1 (simple mean field theory) can give results ap- 
preciably disparate from MC simulation results, larger 
clusters (m = 2,3,4,...) significantly improve conver- 
gence to the correct mean density level. Moreover, sim- 
ple MFT (Eq. [2|) is exact in the symmetric, equilibrium 
limit fc_|_ — > fc_, where the moment equations are closed. 

Now consider a lattice-associated particle that can be 
moved by the nonequilibrium fluctuations inherent in 
the substrate. Simple kinetic rules are defined in Fig. 
[3] Without loss of generality, we also assume that the 
particle does not have intrinsic fluctuations. Since its 
stochastic motion is caused only by the asymmetric do- 
main walls fluctuating past it, we must determine the 
domain wall probability in the frame of the moving par- 
ticle. The thick solid black curves in Fig. [^b show the 
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FIG. 3: The kinetic rules for domain wall-driven particle mo- 
tion. A particle is pushed forward with probability /+ when- 
ever a . . . 100 . . . domain wall tries to move past it from left to 
right. Similarly, the particle moves backward with rate 



mean hydrolysis level, determined by MC simulation, in 
the particle frame. The particle position is arbitrary, but 
the hydrolysis level near it differs significantly from the 
uniform bulk away from the particle (or in the laboratory 
frame). As we follow the stochastically driven particle, 
the mean hydrolysis level (ai) just before (after) it is 
higher (lower); the particle statistically moves ahead of a 
domain wall, spending more time ahead of it. 

The mean velocity and dispersion of the driven particle 
are computed from 



F = fc+Z+Qioo - fc-/-Qooi and 



(4) 



where Qioo and Qooi are the steady-state probabilities 
that the segment of three sites centered about the driven 
particle is in the indicated configuration (cf. Fig. [3]). In 
order to use mean field theory to compute Qioo and Qooi i 
we must use either moment equations or a finite-sized 
mean field transition matrix in the moving, fluctuating 
frame of the transported particle. In analogy to Eq. [2 
we can consider the evolution equation of the first mo- 
ment (cTf) of the hydrolysis level at the site of the driven 
particle. In addition to the state transitions represented 
by rates k± (hydrolysis) and ko (dehydration) , transition 
terms also arise from motion of the driven particle: 



dt 



fc- - (k- + ko){ai) + A{(7e-i) - A((t^o-£_i) 



(5) 

where p± are the effective forward and backward hopping 
rates of the particle. Since the particle moves only via 
motion of domain walls defined by Qioo and Qooi, 



P+ = k+f+QioQ = fc+/+o-f_i(l - <Tf)(l - ae+i) 
P- = k-f-QoQi = fc-/-(l - <7e-i){l - cre)ae+i. 



(6) 



Upon substitution of p± in Eq. [6] into Eq. [5l and 
neglecting all correlations, we obtain a simple, single- 
site, moving-frame mean field approach to finding (a) ~ 
(ct^) « (c£±i) as a root of a cubic equation. This sim- 
ple mean field solution is only in qualitative agreement 
with MC simulations. In analogy to the finite-sized seg- 
ment approach implemented through the transition ma- 
trix M, we can also improve the simple moving-frame 
mean field theory by considering a sliding window of sites 
always centered about the driven particle [lH. Within 
these sites, the configurations are explicitly enumerated, 
transitions involving sliding the segment as it follows the 
driven particle are included, and steady-states are found. 
In our subsequent analyses, we use an m = 5 site seg- 
ment that yields sufficiently accurate results for the pa- 
rameters explored. Henceforth, we rescale time in units 
of kg^ (normalizing all rates with respect to ko), and set 
/± = 1. 
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FIG. 4: (a) The mean velocity V as a. function of the forward 
hydrolysis rate k+ for various fc_. (b) The dispersion D for 
the same values of fc_ . The symbols mark the values obtained 
from MC simulation and the solid curves are results from a 
5— site finite-segment mean field approach. 

Fig. [4^ shows the mean velocity V derived from MC 
simulations and from finite-segment mean field theory 
applied in the moving frame of the convected particle. 
The agreement between MC simulation and the finite- 
segment MFT is quite good provided fc_ -5^ 0.1. For 
k- > fc_|_, the hydrolyzed domains grow backward at a 
higher rate than forward, and the mean velocity < 0. 
The velocity is positive and increases once k-^. increases 
past fc_. However, if k-^ becomes too large, V decreases, 
despite an increase in the hydrolysis asymmetry along 
the track. This behavior can be understood by consider- 
ing Eq. [Hand Fig. [S] When both k± <^ 1, dehydration 
dominates, domains are quickly dissipated, and the par- 
ticle is kicked by a rare 100 domain wall as shown in Fig. 
[5^. Increasing values of k+ > fc_ modestly increases the 
asymmetry and hence V. However, if fc+ is too large, the 
hydrolyzed domains merge into each other, as depicted 
in Fig. [SJ;, reducing the domain wall density Qioo and 
ultimately V. Thus, there is a value of k+ that gives an 
optimally combined domain wall density and nucleation 
asymmetry, as shown in Fig. [Sb, resulting in a maximum 
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FIG. 5: The qualitative regimes that lead to maximum in 
mean velocity Vmax- (a) If fc± <^ 1, most of the lattice re- 
mains unhydrolyzed and there are few domain walls to push 
the particle, (b) When k+ > fc- appreciably, but not too 
large to destroy domain walls, and k- < 1, hydrolysis levels 
are intermediary, and the mean velocity is maximal, (c) For 
k+ or k- extremely large, the hydrolyzed domains coalesce, 
diminishing the number of domain walls. Even though the 
rate k+ is large, the quantity fc+Qioo slowly diminishes. 



mean velocity Vmax- 

For all values of fc_, an extremely large fc_|_ will ul- 
timately decrease the particle velocity. Although A = 
— /s_ might be large, Qioo decreases sufficiently that 
V decreases. The decrease of the mean velocity V in 
the fc+ — > 00 limit can be determined by considering 
a "virial" approximation where only transitions among 
111, 110, 101, and 100 need be considered. Nearly all 
sites are hydrolyzed and we use s « 1 in the mean 
field approximations. Since 100 configurations along the 
lattice are rare and spaced far apart, we can consider 
each triplet of sites as independent and find Qioo ~ 



0(k 



-3^ 

+ 



Therefore, V 



2{ko/k+) 
— > 00. 

Fig. shows that the driven particle dispersion is 
maximal near the maximum \V\. This is a property of 
particles with no intrinsic hopping that move stochasti- 
cally only when they are driven by domain walls. 

In Fig. [S] we explore Vmax and the value of fc-|_ that 
yields Vmax (designated fclj.) as functions of /c_. We find 
Vmax and k"^ by accurately fitting the MC simulations of 
V in Fig. 0] with Padc approximants and indicate these 
optimum values in Fig. [S]with open circles. For large 
the hydrolysis level is high and domain merging prevails 
throughout the lattice (Fig. [5j:), requiring ever-higher k*^_ 
to arrive at a Vmax that is smaller. Lowering to inter- 
mediate values increases both the asymmetry k^—k^ and 
delays the onset of domain merging, giving rise to larger 
Vmax at smaller values of fc^. However, for k- < 0.1, 
the hydrolyzed fraction and the onset of domain merg- 
ing become insensitive to fc_ and the optimal value of k+ 
asymptotically approaches fc^ ~ 3.3. The m = 5 finite- 
segment MFT results are plotted as solid curves and for 
A:_ > 0.1 are in close agreement with those from MC 
simulations. Note that for the case k^ before the 
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FIG. 6: (a) The maximum mean velocity Vmax attainable and 
(b) the value fc^ that yields this maximum velocity for given 
k-. Results derived from MC simulation (circles) and 5-site 
finite-segment MFT (solid curves) compare favorably. The 
dashed curves are numerical fits to the MC simulation data 
given by Eqs. [T] 



system size iV — > cx), cr^ = Vi is an absorbing but un- 
stable condition accessible only by large deviations from 
the stable steady-state that are exponentially unlikely in 
the system size. Nonetheless, for a large system a well 
defined steady-state can be found and Vmax and k^ can 
be accurately fit using 



0.315 



kl 



fc_ + 0.469 

1.24 + 2.67fc_, fc_ > 1, 



(7) 



and Vmax « 1.68exp(-1.84fcO-29)^ /;* ~ 3.3 for fc_ < 0.1. 
These expressions are plotted as dashed curves in Fig. 
inland provide accurate, universal approximations on the 
maximum velocity possible Vmax as a function of fc_ , and 
the value of fci^ required for maximal particle velocity, for 
each fc_. 

In summary, we have presented a paradigm for 
substrate-driven particle motion which has been a rel- 
atively understudied mode of subcellular transport. Our 
model captures the salient aspects of hydrolysis waves 
and exhibits rich transport behavior. Specifically, we find 
short-ranged state correlations, allowing us to accurately 
compute nonequilibrium steady-state particle velocities. 
For fixed backward hydrolysis rate A:_ , the velocities show 
a maximum as a function of the forward hydrolysis rate 
k-i-. The value k'^ approaches its minimum near ~ 3.3 
when k- < 0.2. The maximum velocities and the asso- 
ciated rates can be accurately described by the simple 
universal fitting equations [T] 

Additional details such as slippage (/± < 1) and an 
external load can be readily implemented. A load on the 
particle would bias the motion of the motor backwards 
and impart a negative drift velocity in addition to that 
shown in Fig. |4^. The force-velocity relationship follows 
directly from the functional form of the force-dependent 
drift. 
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